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ABSTRACT 


A review is presented of the results obtained in the research 
program of microwave scattering properties of snow fields. The program 
is a cooperative effort between a group at the University of California, 
Berkeley (D. J. Angelakos, F. D. Clapp, K. K. Mel, and S. Coen) and 
NASA-Ames (W. 1. Llnlor) . A goodly portion of the research was under- 
taken at the Central Sierra Snow Lab and with the active cooperation of 
the U.S. Forest Service (J. L. Smith). 

Experimental results are presented showing backscatter dependence 
on a) frequency (5. 8-8.0 GHz), b) angle of incidence (0-60 degrees), 

c) snow wetness (time of day) , and d) frequency modulation (0-500 MHz) . 
Theoretical studies are being made of the inverse scattering problem 
yielding some preliminary results concerning the determination of the 
dielectric constant of the snow layer. 

The experimental results lead to the following conclusions: a) 

snow layering affects backscatter, b) layer response is significant up 
to 45 degrees of incidence, c) wetness modifies snow layer effects, 

d) frequency modulation "masks" the layer response, and e) for the 
proper choice of probing frequency and for nominal snow depths, it 
appears to be possible to measure the effective dielectric constant and 
the corresponding water content of a snow pack. 


I. INTRODUCTION 


The mapping of ice and snow and soil moisture is an important 
project because such maps provide essential information for water 
resources management and agricultural or forestry yields. Mapping by 
microwave remote sensing has recently been identified as an important 
technology in addition to the optical and Infrared techniques. There 
are two branches of microwave remote sensing techniques currently being 
Investigated for airborne mappings, namely, the radiometry and radar 
remote sensing. The radiometry is a passive sensing technique which 
measures the emissivity of the target subject to illumination by natural 
sources. The radar remote sensing is an active sensing technique which 
measures the echoes of the target subject to illumination by an airborne 
microwave source. There are advantages and disadvantages in both 
approaches, which have been thoroughly discussed in reference [1]. To 
date, the microwave radiometry appears to have more promise for success- 
ful applications than radar remote sensing. That, however, does not 
necessarily mean that radar methods are basically Inferior in remote 
sensing. Most probably, that situation just reflects the amount of 
research effort devoted to these techniques. The art of radiometry had 
been well developed by the radio astronomers before it was brought to 
use in microwave remote sensing; there have been well-documented theories 
relating brightness temperature to the atmospheric composition and 
temperature profiles. On the other hand, radar remote sensing is a 
rather young technology. So far, research effort on radar remote sensing 
is mostly concentrated on instrumentation; it lacks the corresponding 
solid theoretical foundation which supports radiometry. 
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II. ERL*s UNIQUE EXPERIENCE AND QUALIFICATIONS 

IN MICROWAVE REMOTE SENSING OF SNOW 

The Electronics Research Laboratory of the University of California 
at Berkeley has been supported by NASA under contract NASA-Goddard Grant 
NSG-5093 to study active microwave r^ote sensing. The total support 
for the research is $61,000 for two years. With the limited funds and 
time, the ERL team, in cooperation with Dr. William I. Llnlor of the 
NASA-Ames Research Center, was able to obtain rather significant results, 
both theoretically and experimentally, which are basic to radar remote 
sensing. A report of these results is contained in the paper by 
Angelakos, which is attached in the appendix of this proposal. These 
experiments were performed in the frequency band 5.8-12 GHz. An inter- 
esting discovery is the strong dependence of the backscatterlng of the 
snow pack on time of day. The phenomenon was later confirmed by several 
other groups (Curie et al [2], Ulaby et al [3]). Some typical back- 
scattering results are shown in Fig. 1. It is noticed that considerable 
variations were recorded at 15 min. intervals. Due to the time depen- 
dence of radar backscatterlng of the snow pack, empirical correlation 
between measurements and snow properties is difficult to establish. 

That is one of the reasons we feel that the optimum frequency band for 
snow pack remote sensing is yet to be determined. We propose for future 
research to investigate radar backscatterlng starting at lower frequen- 
cies such as 100 MHz and proceeding upwards in frequency. 

On the theoretical front, the ERL team has formulated the integral 
equation describing the electromagnetic scattering from layered media in 
a form in which only the amplitude of the backscatterlng is needed to 
predict the dielectric constants of the layers. The results of a 



OepfAdfftCf Oft Tiao cf Oof 





f a» 


». 0- 

H • -- 

iQddcr 

MofCh ('♦ * T 
( U pm 


Backscattering of a snow pack as a function of 
frequency at 11:30 A.M. , 12:09 A.M. and 1:13 I’.M. 



typical inverse scattering are shown in Fig. 2. It is noted that the 
inverse algorithm we have used predicts the layered media quite well 
even with only three (3) measur^ents of the backscattered amplitude. 

The details of our formulation and algorithm will be presented in the 
following sections. 

III. GENERAL BACKGROUND 

(a) Measurements and Instrumentation 
Measurements of microwave backscatterlvig by snow pack have been 
made by many investigators. However, with no exception, all previous 
efforts were channeled towards finding the electrical properties of the 
snow, assuming the snow pack is one uniform layer. The objective of 
ERL's snow measurements has been to understand the field environment 
and to develop instrumentation for untlmately collecting data for an 
Inverse scattering theory which predicts the profile of the snow property. 
Toward that end, it was believed that the repeatability of the measure- 
ments should be of the greatest Importance. It was to our great 
surprise that the backscattering of snow pack is not repeatable, even 
at an interval as short as fifteen minutes. The experiments were made 
at the Central Sierra Snow Laboratory at Tahoe, performed by a te^un of 
the ERL, in cooperation with Dr. William Linlor of NASA-Ames and the 
U.S. Forest Service (J. L. Smith). The instruments involve horn 
antennas mounted on towers as shown in the picture in Fig. 3. The 
results obtained at 15 min. Intervals at 8.2 GHz and 10 GHz are shown 
in Figs. 4 and 5. It is observed that the variation of the backscatter- 
ing of the same snow pack is so unpredictable that the simplistic 
approach of predicting snow electrical property by comparing its radar 
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backscatterings with a set of known signatures is hopelessly an 
Inadequate model. 

With the above understanding of the dynamic property of snow pack, 
we feel that two things should be done in order to affirm the feasibil- 
ity of radar remote sensing of snow packs, namely, first, further 
measurements be done at much lower frequencies, such as 100 MHz, or 
even as low as 25 MHz, and secondly, concrete evidence of the applica- 
bility of the Inverse scattering technique be established. 

The purpose of lowering the frequency of observation is two-fold: 
to determine whether the time dependence of the backscatterlngs from 
snow is a surface effect, and to condition the instrumentation for 
collecting data towards the use of an Inverse scattering algorithm. If 
the time variation of a snow backseat ter Ing was a surface effect, lower- 
ing the freque.icy should reduce the variation greatly. Furthermore, 
for successful prediction of the snow properties by inverse scattering, 
the signal must penetrate the entire anow pack. It would be impossible 
to predict the snow properties for regions which the field cannot reach. 
In many field situations, the snow could be over 20 feet deep, hence it 
woulf’ be unlikely for a 1 GHz wave (X“30 cm) to penetrate a wet or packed 
snow layer of that thickness. 

It should be eaphaslzed here that the time variation of radar 
backscattering Is not a deterrent to the inverse scattering approach, 
because the backscattering dynamics are most probably due to the effect 
of changing snow profile, l.e., the backscattering data are fairly 
sensitive to the snow prorile constituency. Since the end result of 
the inverse scattering technique is to find the snow profile, the 
sensitivity of the measuranents indeed offers prospects for better 


resolution. We shall discuss the backgroimd for Inverse scattering In 
the following section. 

(b) Theory of Inverse Scattering 
(1) Introduction 

Inverse scattering pr> blems in science and engineering concern 
the detemln<itlon of the properties of some Inaccessible region from 
observations away from the region. The observations are often taken as 
the response of the Inaccessible region to some source, which may be 
natural or artificial. It may be located within or on the boundary of 
the inaccessible region. 

The Inverse scattering theory falls Into two distinct parts. One 
deals with the Ideal case In which the observations are assumed to be 
known exactly and as dense as required, whereas the other treats the 
practical problems that are created by Incomplete and imprecise data. 

To apply these Ideas, there must be a valid mathematical model that 
describes the physical behavior of the Inaccessible region. The model 
is a plane stratified half-space whose properties are characterized by 
an inhomogeneous, linear and Isotropic dielectric. 

The equation which describes the model is a one-dimensional wave 
equation. 


+ k^ e(x)l E(k,x) - 0 
^dx ^ 


—O' < X < 


( 1 ) 


where e(x) is the dielectric profile and E(k,x) Is the electric field. 
We have assumed the wave to be polarized parallel to the half-space, 
hence the polarization direction of the electric field Is omitted. 

The conventional interest in eq. (1) is to oolve the electric 
field everywhere, given the Incident field and the dielectric profile 
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e(x), which, In the remote sensing termlnolo^ r is known as the forward 
problem. There are many different ways the forward problem can be solved 
economically and accurately. In the ensuing discussion, we shall assume 
that all forward solutions of (1) are available for any given t.(x). 

The inverse problem of (1) can be posed as the following: 

Given the Incident field E^(k,x) and the reflected field E^(k,d) at 
a fixed position d, f id the dielectric proixie of t(x). 

The inverse problem wheie E^(k,d) is known for all k, (0 < k < 
can be solved by Gelfand and Levitan theory [4), which also proves the 
uniqueness of the inverse scattering problem. However, Gelfand and 
Levitan theory is limited to real l(x). 

If the medium is lossy, i.e., complex c(x), and the information on 
the reflected field is Incomplete, i.e., E^(k,d) is known at a finite 
number of frequencies, the Backus and Gilbert's iterative technique can 
be used [ 6] . We shall discuss the formulation and application of the 
Backus and Gilbert method in the following paragraph. 

(2) Iterative Procedure 

Consider two different wave e<; iatlons of the type of eq . (1), 


— r + k c (x) 
j 2 n 

dx 


E^(k,x) 


n = 1,2 


( 2 ) 


The scattering geometries are shown in Fig. 6. Multiplying the equation 
for n “ 1 jy E^, and that for n = 2 by subtracting the resulting 
en-iMtions and integrating over the snow region, we get 

oo(k) ■ yi’ \ >5t(x) Ej^(k,x) E^Ck.x) dx (3) 
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where 


5p * and 6c = ~ 


Let £j^(x) be a given profile, then by solving the direct problem 
numerically, E^(k,x) and Pj^(k) are known functions. If p 2 (k) were 
given from measurements, then eq. (3) is a nonlinear equation for 
£ 2 (x). One way of solving the nonlinear equation is by developing an 
iterative scheme, based on Newton’s like method in functional space. 

The linearized equation to be solved at each step of the iteration scheme 
is given by 


6p(k) 


_k_ 

2j 


5e(x) 


E (k,x) Jx 
1 


(A) 


k 2 

where -ry E (k,x) is the Frechet derivative. 

1 

Normally, eq. (4) requires measurements of p 2 (k,d) for both 
magnitude and phase. Phase measurement is, in general, difficult to 
Implement. To circumvent the requirement for phase measurement, we have 
reformulated eq. (4) by a transformation to. 


g^ = 1 f(x) L^(x) dx i = 1,2,...,N (5) 


Ipjdtj) 

where = log 

r ^1 

and L^(r) = 


, f(x) = - e^(x) 


(k ,x) 
1 ^ 


where the wave number k is 


discretized and the medium is assumed to be lossless. N represents the 
total number of measurements. Note eq. (5) requires only the magnitudes 
of backscattered signal. 

(3) Resolving Power Theory 


Suppose that after finite iterations a profile has been determined 


which fits the data. The question is then, how good is the resultant 
profile? This question can be answered by Backus and Gilbert theory 
[6], as follows. 

After the iteration of (5) converges, which means g^^ = 0, for 
i “ 1,2,...,N, we have 


f ^ 

= \ E]^(x) L^(x) dx » 

Jo 


h 

e^(x) L^(x) dx 
o 


i « 1,2 N. (6) 


where and L^(x) are known. Multiply the equations by a^(x) and 5jum 
over i to obtain 



< £ 2 ^ “V ^2^*^ A(z,x) dx 

Z Jq 


(7) 


where 


N 




-Z 

i=l 


“i ®i^^J 


( 8 ) 


N 

A(z,x) = Xw a (z) L,(x) 
i«l 


(9) 


If the functions selected that A(z,x) resembles the Dirac 

Delta function 6(z-x), then <£ 2 ^ may be viewed as a local average of 

z 

at z. Further, the local average < e-> is a common property of all 

^ z 

the solutions that satisfy the data. The functions a^(z) are determined 


by 


min 




(z-x)^ A^(z,x)dx 


( 10 ) 


o 

.h 


i.t. ^ A(z,x)dx = 1 
■^o 


Upon introducing a Lagrange multiplier, the minimization problem is 
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reduced to a set of N+1 linear algebraic equations. Once the coeffi- 


cients a^(z) are determined, the local average ( is given by 

z 

eq. (8) and the averaging kernel is computed by eq. (9). 

Another such "Deltaness" criteria is the Dirichlet Kernel. 

min. 1 [ 6(r,-x) - A(z,x)]^ dx 

*^o 


( 11 ) 


s.t. 1 'V(z,x)dx * 1 
"'o 

Fig. 7 shows the averaging and the Dirichlet kernel. It is clear 
that there is a trade-off between resolution and side lobes. 

In case the data is erroneous, one minimizes the spread of the 
kernel and the error in the solution. This leads to a trade-off between 
resolution and certainty. 

(4) Results with Artificial Data 

Fig. 8 presents results for a typical snow profile; the results 
sliown in Fig. 2 have been obtained with data at 3 frequencies only, 
whereas the results shown in Fig. 8 have been obtained with data at 7 
frequencies. The approximate results are in good agreement with the 
exact results for smooth profiles, as is evident from Fig. 2. The poor 
resolution near the edges of the ice layer (Fig. 8) may be improved by 
using more data. It will be noted that the computation time was about 
7 seconds on a CDC 7600 computer for the profile reconstruction shown 
in Fig. 8. (computation cost “ 2.5 $). 

(5) Best Solution and Error Estimates 

The linear Fredholm eq. (5) does not have a unique solution. 
Further, the data g may have some noise. Miller and Viano [5] have 
shown how to obtain solutions to eq. (5) with prescribed bounds and 
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beat possible error estimates. They suggest to find f which satisfies 


iLf-gl _< a 

( 12 ) 

Ifl < $ 

where a Is a small number and 0 Is a fixed number. Obviously a Is 
related to the noise in the data and 0 may be obtained from the minimum 
and maximum values of the dielectric constant of the medium, which are 
assumed to be known. They further show that f can be obtained by 
minimizing the functional 

J(f) - iLf-gl^ + (|)^ Ifl^ (13) 

which leads to 

Cf - [L*L-^(^)^ I] f - L*g (14) 

P 

* ot 2 

Eq. (14) tan be inverted since L L is symmetric and (^) > 0. The solu- 

P 

tlon of (14) is obtained by expanding f in terms of the eigenfunctions 

* f 1 

of the syna&etrlc operator L L. As suggested In ref. 151, the error 
bound can be extracted from C ^ in (14) . 

IV. FUTURE RESEARCH 

Additional research needs to be conducted, both experimental and 
theoretical. Both are discussed In the following. 

(1) Instrumentation and Measurements 

(a) Backscattering at Low frequencies 
The type of measurements we have been doing at the snow laboratory 
at Tahoe should be continued and extended to lower frequencies. Also, 
the results should be recorded digitally for direct application In 
inverse scattering. 
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For operation at lower frequenclea down to about 100 MHz the 
physical arrangement must be changed. Considering the wide percentage 
bandwidth needed, log periodic antennas could be us^. For a range of 
100 - 1000 MHz, the dimensions of the log periodic antenna will be on 
the order of 1-1/2 meters wide atUl 2-1/4 meters long, which would be too 
large to be mounted on a frame. It is believed that satisfactory opera- 
tion can be obtained by suspending the antennas from a cable between 
two poles or trees with guy stabilizers. 

Due to cyclic variation of gain with frequency and the added cable 
loss at the high frequency end of the 10 ;1 b<uid, it is recommended that 
a storage-gain notmallzer such as the HP 9750A be used. The response 
curve of the system as a function of frequency can be generated and 
stored in this Instrument. These data can then be used to automatically 
normalize the measuranents , eliminating the system error. 

(b) Measurements in Controlled Environment 

Measurements in controlled environment are a necessity in all 

experimentation. It is particularly true in the snow measurement 
because at the snow laboratory snow is only available in the winter 
season, and the instruments muiit be calibrated and ready when the snow 
season arrives. A sand pit can be used to simulate layered media, using 
polyethylaie sheets to separate layers of different wetness. The 
experiments should Include the following measurements: 

(1) Calibration of instruments to be used in snow measurements . 

(2) Gathering data for Inverse scattering . This involves layered 
wet sand. 

(3) Studying the effect of rough surfaces . 

(c) Ground Truth Collection 
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(I) latftrferoaeter for dielectric conatMit and 
attaauatlon meaaurementa 

One convenient method of obtaining anow truth for comparlaon with 
the varloua meaaurem«it8 has been the use of an X>band (9.3 Hs) Inter- 
ferometer. This set-up can give the dielectric constant and attenuation 
directly at Its frequoacy of operation. The major problem oil operating 
such equifsaent at VHF frequencies Is the antenna » which illuminates the 
staple. However » It appears that comer antennas can be designed to 
fulfill this function in the range 500 - 1000 MHz. 

(II) Attenuation In snow 

While attenuation Is fairly easily measured in the Interferometer 
by sampling, measuranents in situ have not been so simple, due mainly to 
layering effects. A preliminary test last year made by burying a 
receiver in the snow and illuminating it vertically frena the top appears 
to be a more satisfactory procedure. For this type o£ measurement, the 
receiver antenna and detector are best placed on groimd level in a 
weatherproof enclosure before the snow falls. As the snow builds up 
during the winter, a transmitting unit can be positioned above as 
required. Thus no disturbance of the snow need occur. 

(2) Theory 

(a) Algorithm for Lossy Media Inverse Scattering 

The Inverse scattering results shown In Fig. 8 are for media of 
real dielectric constants. Future research will extend whe inverse 
scattering technique to complex media. For this case, eq. (5) will be 
replaced by 

*1 " 3 L^(x)dx + \ s(x) T^(x)dx 1 - 1,2 N (15) 
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«^ere L^(x) and T^(x) are known and f(x) and s(x) are sought; here, 
s(x) ■ 02 (x) - o^Cx) and o(x) denotes the conductivity profile. 

The resolving power of the data in this case may be obtained by 
constructing averaging kernels which resemble 5^^ 6(x-^), where is 
the Kronecker delta and fi(x-C) is the Dirac delta function. Such 
kernels can be constructed by using the criteria suggested by Backus and 
Gilbert [ 7J . 

(b) Error Estimate 

The best error estimates of the complex medium Inverse scattering 
can also be achieved using the Miller and Vlano (5) method, except now 
we have to replace (12) by 

iLf + Ts - g _< a 

Ifl < e (16) 

Isl ^ Y 

(c) Inverse Processing Experin«ital Data 

The ultimate goal should be to process the data collected in the 
field and c(X&pare the predicted snow profile with ground truth. 

V. CONCLUSION 

The remote sensing of a snow pack Is a ccxaplex technical problem. 

So far, neither the passive microwave technique (radiometry) nor the 
active microwave technique (radar) can claim to be able to solve the 
probl^ completely. Our effort has been unique in that its objective is 
to join experimentation and theory of Inverse scattering. Our future 
goal is to predict the electrical properties of snow pack by remote 
sensing, which we believe Is attainable. However, the useful Information 
of the snow is its water content, (both liquid phase water and water 
equivalent) which is not just a function of the electrical property. 
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Ttt^peratuttt, for exuiple, is an important parameter, affecting the 
relation between the electrical properties of snow and its water 
content. It is possible that the information gathered by the radlometry 
method will provide the temperature profile and together we have a good 
chance to eventually predict the water content of sncn^. We are confl- 
d«it that the fruit of our future research will bring the remote 
spaing technology to a new level of sophistication and push the 
frontier much closer to the final goal. 

The purpose of this investigation is to provide radar rmmite 
sensing with the basic theoretical and experimental supports \dilch it 
so badly needs. The objectives of future research can be susBaarlzed as 
follows: 

(1) To idratlfy the frequency bands most advantag<»}U8 for remote 
sensing of ice, snow fuid soil sujisture. 

(2) To establish the thet^retlcal relation between the radar 
return signals (at selected frequencies) and the dielectric 
profiles of the target areas. 

(3) To design experiments in controlled environmrat to test the 
theory, and the statistical properties of the target area. 

(4) To design experiments for measurement in field conditions and 
to develop techniques of gathering reliable ground truth 
information. 
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APPENDIX 


MICf^OWAVE SCATTERING PROPERTIES OF SNOW FIELDS 
D. J. Angelakos 

Electronics Research Laboratory 
University of California 
Berkeley, CA 94720 

ABSTRACT 

A review is pr^ented of the results obtained in the research program of 
i^crfwsve scattering properties of snow fields. The program is a cooperative effort 
bet.w’^n a group at die University of California, Berkeley (D. J. Angelakos, F. D. 
Cl« 4 »y, K. K. Mei, and S. Coen) sid NAS A- Ames (H. I. Llnlor). A goodly portion of 
the research was undertaken at the Central Sierra Sm>w Lab arKl with the active 
co<^eratioa of the U.S. Forest Service (J. L. Smith). 

Experimental results aie presented showing backscatter dei>endence on a) 
frequency (5. 8-8.0 (^r), b) angle of incidence (0-60 degrees), c) snow wetness 
(time of day) , and d) frequency modulation (0-500 M*4z) . Theoretical studies are 
being s^e of the inverse scattering problem yielding some preliminary results 
ccmceming the determination of the dielectric constant of the snow layer. 

The experimental results lead to the following conclusions: a) Snow layering 

affects backscatter, b) L^er response is significant up to 43 degrees of Incidence, 
c) Wetness ^>difles snow layer effects, and d) Frequency modulation "masks" the 
l^er response. 


Introduction 

At the Electronics Research I.aboratory, University of California, Berkeley, the 
research program of microwave scattering properties of snow fields along tvo 
directions. Oi^ Is the experimental determination of the microwave properties of 
snow fields and the other is a theoretical program of mathematical modeling ami the 
Inverse scattering probl^. Hie ultimate goal of the project is the coalescence of 
the two into a useful diagnostic technique for the detexminatlon of water content 
in anew. 

E xperimental Portion of the Project - D. J. Angelakos, F. D. Clapp, W. I. I ini or 
a»l J. L. Smith: The primary objective of the experimental program was to determine 

the effect of snow couditlons on the asount of backscattered elec romagnetlc energy. 
In addltl(xi the investigation was to seek out any effects on th.' i> kscatter due to 
SfUMf layering and to determine the dependence on frequency and frequci.cy deviation. 

In proceeding towards the objective several additional phenomenological effects 
appeared *dilch added insight to backscatter signal response. What follows is a 
saaqillng of results idiich liulicates the partial attainment of cur goal. 
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EquitOTent ; Most of the data were taken for swept-frequency ranging from 5.8 to 
8.0 GIz. The sweepers were the "Alfred" and the "Hewlett-Packard" sources. These 
were connected with coax cable to horns having nominally 20-dB gain. The receiver 
horns were coupled to crystal detectors and HP SVTR meters, whose cxitput was plotted 
on XY plotters, on a logarithmic scale. For some measurements a converter unit was 
employed (Wlltron manufacturer) that permitted a uniform dB scale for the Y-axis. 

For the 13.5 GHz measurements a fixed- frequency solid state oscillator was 
used (Varian manufacturer) . Time was not available to take measurements at other 
frequencies, either flxa) or swept. 

Sites I Two sites were used: one is called the "tower site" which is on the roof of 

the snow hut, for which only two fixed angles (39“ and 45“) could be used; the other 
is called the "ladder site" where all desired angles could be obtained. 

The limited available time — dictated by the snow depth — for measurements 
and also the limited manpower made a choice necessary between taking data or ground 
truth (i.e., snow electrical characteristics). Under the circumstances the ground 
truth had to be limited to: 

a. Density profile of the snow 

b. Visual inspection of layers in the snow and granularity 

c. Dielectric constant measurement of a few samples, near the top and middle 
of the snowpack 

d. Wetness measurement of th^ , > samples 

Discussion of Experimental Results : The solid line of Fig. 1 shows the backs catter 

in dB versus frequency in the range of 5.8 to 8.'/ GHz, at the fixed angle of 39® 
from the vertical, and a snow depth of 70 cm. The electric field vector of the 
horns (20 dB gain each for transmitter and receiver) was parallel to Lhe earth 
surface. 

The various peaks are produced by layering effects in the snow. The peak 
values are about 25 dB greater than the minima. 

To then determine vhat effect frequency devlatiot’ (more specifically, frequency 
modulation) had on the measurement the signal source was set to the frequency for 
each maximum and minimum, but a modulation of +100 MHz '^aj superimposed for the + 
points; and a modulation of +200 MHz was superimposed for the 0 points. Evidently 
the effect of frequency modulation is merely the averaging of the unmodulated 
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response over the bandwidth represoited by the modulation. 

In the course of taking measur^ents , %diether at normal Incidence or otherwise 
It fNis observed that the backs catter Ing data depended on the observation tl«"*» aid 
temperature . Figures 2» 3 and 4 show the backscatter response at normal Incidence. 
The experimental condltioi^ for the three sets of data were identical except for the 
time taken and accompanying aid>lent temperature. Variation of backscatter with snow 
conditions* at 13.5 Qlz and Incidence angle of 39” Is shown In Fig. 5. 

The decrease In backscatter shortly before noon was Initially observed by us In 
Harch 1976 for a snowpack 150 cm deep* and has been confirmed at other frequencies 
by several grotq>s (Currie et al., 1976* Ulaby et al.* 1977). Evidently the Increase 
In surface snow wetness by solar heating Is the reason for the backscatter decrease* 
totaling about 15 dB. The air temperature variation Included is useful In 
correlation the phenomenological effects. 

Snow wetness produced during a few hcnirs of sunlight can significantly affect 
the frequency response. March 5* 1977, when the meeisurements were taken* was a 
cl«ar, warm day with a taiq>erature of about 10” C between noon and 16:00* so the 
snow was essentially saturated with water In the upper surface. 

Figure 6 shows typical data demonstrating that variation In backscatter as a 
function of frequency can be observed even for other than normal Incidence (39”) . 
Significant backscatter data has been obtained for incidence angles approaching 45”. 

A continuous variation of incidence angle from the vertical is shown in Figure 
7* giving the backscatter in relative dB for the fixed frequency of 11.8 GHz (no 
modulation* with snow depth of 53 cm. The backscatter Is down about 20 dB at the 
angle of 40” from the vertical. 

T heoretical Portion of the Project - K. K. Mel and S. Coen: It Is suggested that 

Certain useful properties of the snow-pack may be approximately determined by using 
inverse scattering theory; the Inverse problem is that of determining the properties 
of inaccessible snow pack from the back-scattered electric field due to normally 
incident plane waves at varies discrete frequencies. 

Let the back-scattered frequency response of the snow-pack be taken as the 
reflection coefficient amplitude |p(wj^)|, 1 “ 1,2,...,N. Neglecting for the moment 
the experimental errors* the question Is: how can we learn as much as possible 

about the characteristics of the snow-pack from a set of numbers |p(u)^)|, 
i • 1*2,..., N? Alternatively what do the numbers [p(u)^)|, i = 1,2,...,N tell us 
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about the dielectric constant and coiKluctlvity of the snov>pack? The density and 
wetness of the snow may be then approximately determined by using Weiner's theory of 
dielectric mixtures or other available experimental relations. 

In what follows, we show that there is enough information in the numbers 
1 p(w^)|, 1 “ 1,2,...,N, to construct local averages of the snow profile. 

Here the experimental errors are neglected and a mathematical model that 
describes the physical behavior of the snow-pack is chosen. The simplest model 
that we choose, is that of a plane stratified snow-pack, whose properties are 
characterized by a linear, isotropic and nonhoraogeneous dielectric medium, as shown 
in Fig. 8a. Maxwell's equations now show that, for time harmonic dependence 
the electric field E, satisfies the one dimensional wave equation 

3 ^ 2 

+ k q(x,k)> E(x,k) » 0 (1) 

3x 

where k * a)(p e is the free space wave number and q(x,k) * c(x) - j is the 

o 

profile function. Here e(x) is the dielectric constant of the snow and o(x) its 
conductivity. 

Consider two media as shown in Fig. 8b. c. The electric field in the respective 
medium satisfies 

3 ^ 2 

+ kV(x,k)} E (x,k) = 0 n = 1," (2) 

3x 

and the scattered fields are outgoing. Equation (2) can be reformulated into an 
integral equation form tuitable for iterative solutions, 

6p(k)“-^J 6q(x,k) E^(x,k)dx + a(k) 


where we have used the relation = Ej^+6E and 6p = P 2 ~P]^ 
measurement without error. It may be shown that the term a(h) goes 
#6q8^, where 11*1 is the norm. Thus, if 6q is small enough, this 
dropped, and we obtain a linear approximation 


were given from 
to zero as 
term may be 


6p(k) 



E^(x,k) dx 


(4) 
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which is a linear Fredholm int^ral equation, for the determination of 6q from the 
given functions 6p and Once 6q is found, the profile q^ is given by 

However, equation (4) is not exact, so the profile q 2 will not in 
general, satisfy the equation (2). To overcome this, we solve the equation by 
iterations, until conv ergence occurs. In practice few Iter ations will suffi ce. 

In practice, the amplitude of the reflection coefficient is much to be 
preferred than its pha.se, and equation (4) can be rewritten into the form 


®i 


h 

6e(x) Fj^(x)dx 

0 


1 - 1,2,. ..,N 


(5) 


where * log 




Pl(ki) 


5e * e«-e. 


2-h ElOc^.x)}; here 

the freqtiency k is discretized and the medium is assumed lossless. The approach to 
be discussed is however applicable to lossy medium as well. 

We now show how to solve Eq. (5) by the Backus and Gilbert approach. Multiply 
the equation by and sum over i; this results in 

h 


< 6e ), 


io 


6e(x) A(C,x)dx 


where 


( 6e ). 


i*=l 


( 6 ) 


(7) 


is the local average of 6e, and 
N 

A(5,X) » £ aAO fax) 
i»l 


( 8 ) 


is the averaging kernel. 

The coefficients i “ 1,2,...,N are determined, such that the averaging 

kernel resembles the Dirac-Delta function most closely. One way of doing this is by 
minimizing the quadratic form 

( 5 -x)V(C,x)dx (9) 

Jo 
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subject to 


I A(^,x)dx - 1 (10) 

*0 

for each this minimization problem may be reduced to the linear equation upon 
Introducing a Lagrange multiplier. 

Once the coefficients ate determined, the local average < is 

computed via equation (7). 

In Figs. 9 and 10 we present results for 3 different snow profiles; the results 
shown in Figs 9a and 9b have been obtained with data at 3 frequencies only, 
whereas, the results shown in Fig. 10 have been obtained with data at 7 frequencies. 
The approximate results obtained from our analysis are in good agreement with the 
exact results for smooth profiles, as is evident from Figs. 9a and 9b. The poor 
resolution near the edges of the ice layer (Fig. 10) may be improved by using more 
data. It will be noted that the computation cost is about 5 dollars per profile 
on CDC 7600. 

We are currently investigating the influence of experimental errors on the 
resulting snow profiles. 

We find that there is a tradeoff between resolution and uncertainty, which is 
in agreement with the Backus and Gilbert theory. 
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